Double Rosensweig instability in a ferrofluid sandwich structure 
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We consider a horizontal ferrofluid layer sandwiched between two layers of immiscible non- 
magnetic fluids. In a sufficiently strong vertical magnetic field the flat interfaces between magnetic 
and non-magnetic fluids become unstable to the formation of peaks. We theoretically investigate the 
interplay between these two instabilities for different combinations of the parameters of the fluids 
and analyze the evolving interfacial patterns. We also estimate the critical magnetic field strength 
at which thin layers disintegrate into an ordered array of individual drops. 

PACS numbers: 75.50.Mm,47.20.Ma 



I. INTRODUCTION 

Ferrofluids are colloidal suspensions of nano-size ferro- 
magnetic grains in a carrier liquid like water or oil [|. 
The dipole-dipolc interaction between the ferromagnetic 
particles is for moderate volume concentrations rather 
small and ferrofluids hence behave magnetically as su- 
pcrparamagncts. Accordingly in the absence of an ex- 
ternal magnetic field the magnetization of the fluid is 
zero. If a field is switched on the magnetic moments of 
the particles orient themselves along the field direction 
giving rise to a macroscopic magnetization. The notion 
superparamagncts refers to the unusually high value of 
the magnetic susceptibility, \ = 1.— 50, to be compared 
with x — 10 -4 for atomic paramagncts. Hydrodynami- 
cally diluted ferrofluids behave like ordinary Newtonian 
liquids with additional contributions in the bulk and sur- 
face force densities stemming from the interaction with 
the magnetic field. 

Due to the unique interplay between hydrodynamic 
and magnetic degrees of freedom ferrofluids show a va- 
riety of instabilities and pattern formation processes. 
Among the most striking phenomena in this respect is 
the so-called Rosensweig instability in which the flat free 
surface of a ferrofluid becomes unstable when subjected 
to a sufficiently strong vertical magnetic field Q- Al- 
though both gravity and surface tension favor a flat sur- 
face the decrease in magnetic energy for a periodic array 
of peaks and troughs can be large enough to overcompen- 
sate the increase in potential and surface energy. Both 
the linear instability and the details of the pattern for- 
mation as revealed by a weakly non-linear analysis have 
been thoroughly studied @, H H S H ■ 

In the present paper we investigate a sandwich struc- 
ture in which a ferrofluid layer of given thickness is placed 
between two immiscible non-magnetic liquids. The sys- 
tem is prepared such that in the absence of a magnetic 
field the layering is stable, i.e. the lower layers have 
larger densities than the upper ones in order to prevent 
the Rayleigh- Taylor instability. Applying a homogeneous 



external magnetic field perpendicular to the undisturbed 
interfaces gives rise to Rosensweig instabilities at both 
the lower and the upper interface of the ferrofluid layer. 
Due to the non-local character of the magnetic field en- 
ergy these instabilities are coupled with each other. We 
first study the interplay and competition between these 
instabilities within the framework of the linear stability 
analysis. Depending on the parameters of the system one 
interface dominates and "slaves" the other one to its un- 
stable wavenumber or both interfaces become unstable 
at rather similar values of the magnetic field giving rise 
to a competition between the corresponding wavenum- 
bers. This is similar to what occurs in Rayleigh-Bcnard- 
Marangoni convection in systems of two superimposed 
fluids which are coupled viscously and thermally at their 
common interface 
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In order to characterize the patterns evolving from the 
instability we perturbatively probe into the weakly non- 
linear regime by expanding the free energy of the sys- 
tem in powers of the amplitude of the surface deflec- 
tions generalizing the methods developed in 0, 0, IT^| . 
When the amplitude of the surface deformations becomes 
comparable to the thickness of the ferrofluid layer itself 
the layer may be decomposed into disconnected parts. 
Within our non-linear analysis we are able to estimate 
the field strength necessary for such a disintegration to 
occur. Finally by using experimentally relevant values 
for the parameters we point out interesting experimental 
realizations of our system. 



The main difference between our sandwich system and 
the somewhat related problem of a ferrofluid film inves- 
tigated in 01 is the thickness of the ferrofluid layer. For 
a film this thickness is by definition much smaller than 
the wavelen gth of the unstable mode. In the experiments 
reported in |l3[ the film thickness varied between 5 and 
60 fj.ni. The hydrodynamics of the film can then be very 
well described within the lubrication approximation. In 
our system the thickness of the ferrofluid layer is compa- 
rable to the unstable wavelength which is of the order of 
centimeters and correspondingly the full hydrodynamic 
equations have to be solved to describe its dynamics. 
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FIG. 1: Schematic two dimensional plot of a ferrofluid layer of 
depth d with infinite horizontal extension sandwiched between 
two non-magnetic liquids of infinite depth. 



II. BASIC EQUATIONS 

We consider a horizontally unbounded ferrofluid layer 
of thickness d and density p^ sandwiched between two 
immiscible, non-magnetic liquids with densities and 
p^ . The interfaces between the layers are parametrized 
by the functions z = ^ d \x) and z = ^°\x) where for 
simplicity we will only consider one-dimensional interface 
modulations (see fig [I}. It has recently been clarified that 
this situation can be realized experimentally by using an 
oblique magnetic field |l4j . The interface tensions at the 
two interfaces are denoted by and . 

The hydrodynamics of the system is quite generally 
described by the Navier-Stokes equation and the conti- 
nuity equation. However, since the situation of interest 
is a static one these equations can be replaced by the 
pressure equilibrium at the two interfaces. This in turn 
is equivalent to the minimum condition for the total en- 
ergy functional. 

The energy per area in the x-y plane comprises three 
parts, Eh, E s , and E m denoting the hydrostatic, the 
interfacial, and the magnetic energies respectively. The 
first two parts are given by the well-known expressions 

E h = glp { - 1) J dzz + pW J dzz + pW J dzz\ (1) 

^ -oo C (Q) (^) (( d >(x) ' 



E s = (a® yf\ + (d x (W (x)r + <y {d) yfi + (d x (W ( x ))< 



(2) 



where g is the acceleration due to gravity. The brackets 
(...) denote the spatial average along the cc-direction 
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(3) 



The volume density of magnetic energy is of the general 
form [l3| 



-Pa 



Ho 



dH' ■ M(H') 



(4) 



where M denotes the magnetization, Ho the magnetic 
field in the absence of any permeable material, and po is 
the permeability of free space. Assuming a linear magne- 
tization law M = %H of the ferrofluid with the suscepti- 
bility x characterizing its magnetic properties we hence 
find in the present case for the magnetic energy per unit 
area 



C M (x) 



E n 



Mo X 



J dz H(.t, z) ■ H c 



(5) 
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Here H cxt denotes the homogeneous external magnetic 
field produced by the experimental setup and in the ab- 
sence of the ferrofluid and H(x, z) is the actual magnetic 
field in the ferrofluid. 

Subtracting an irrelevant constant the complete energy 
functional of the system can hence be written as 



C (0) W,c (d) (x) 



C (d, (*) 
Mo X 



(p (D _ p (2) )c (0) 2 + ( p (2) _ p(3) ) C (d)=j -f^jdz H(x, z) ■ H cxt + a® yJl + (d x C«»(x))* + a™ yjl + (d x (W (x)f 
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(6) 



The magnetic field has to obey the magnetostatic Maxwell equations 

V • B = and V x H = (7) 
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with B = /io(l + x)H. These equations are completed 
by the following boundary conditions 



lim H(x, z) = H ext e z 

Z — >-±00 



and 



(BO) -B^j.nW 
(HO) -H< 2 >) xn» 
'(B< 2 > -B^-n^" 
(H^-HWlxnl ' 



z=C (d) 



= 



(8) 



(9) 



= £(<>) 



= 



(10) 



0. 



where n^ ' and n^ d ' denote the normal vectors on the 
lower and upper interface respectively. Note that the 
last four boundary conditions have to be fulfilled at the 
free interfaces of the ferrofiuid layer. They hence describe 
the feedback of the interface modulations on the magnetic 
field. Note also that therefore the energy © depends in 
a complicated non-local way on the surface deflections 
C (0) (x) and (W(x). 

It is useful to introduce for each of the three liquid 
layers a scalar magnetic potential and $' 3 ' re- 

spectively. The potentials arc related to the correspond- 
ing magnetic fields by 

H« = -V$ w (11) 

and as a consequence of (Q they obey the Laplace equa- 
tions 



A$ w = 0. 



(12) 



The boundary conditions (JHJ) and (|lt)|) for H and B trans- 
late in the well-known way into conditions for the con- 
tinuity of the potentials themselves and jumps of their 
normal derivatives |l7j . 

It is furthermore convenient to measure all distances 
in units of the inverse critical wavenumber 



1 



r (d) 



( p (2) _ p (3)) g 



(13) 



of the Rosensweig instability on an infinitely deep fer- 
rofiuid layer, all magnetic fields in units of the corre- 
sponding critical Rosensweig field 



(1+X) (2 + x)^(p^-p^)ga(") 



X 2 Po 



-, (14) 



and energies per area in units of a^ d ' . Moreover we in- 
troduce the parameter ratios 



Pi 
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(i) 

(2) 
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a (0) 



P3 



,(3) 

12) 



(15) 
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with i] now characterizing the magnetic properties of the 
ferrofiuid. After rescaling the magnetic potentials ac- 
cording to 



_ ( 2 + __>$(!) 



(1 + X )(2 + X ) ^ 2) ^ (2) 

( 2 + X) ^(3) $ (3) 

xHcxt 

the energy © assumes the dimcnsionlcss form 



(16) 
(17) 
(18) 
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The boundary conditions I© and 11UI) translate into 



^($0) _ $(a))0 x £(«o _ a z ($< 3 ) - $< 2 ') 



1 + ?? $(3) _ $(2) 
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z=C (d) 





(20) 



and 



^($(2) _ $(i))^C (0) - 9 Z ($ (2) - $ (1) ) 



$(2 



l-r, 



z=C<°) 



= 





(21) 



respectively whereas the asymptotic boundary conditions 
(JHJ) acquire the form 

lim d z ^ 3) (x,z) = lim (x, z) = - . (22) 



III. LINEAR STABILITY ANALYSIS 

In this section we study the linear stability of the ref- 
erence state with flat interfaces C,^ = and = d. 
To this end we use the ansatzes 



C (0) = A 1 cos(kx) 
C (d) =d + B 1 cos(fcx) 



(23) 
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for the interface profiles. The corresponding forms of the 
magnetic potentials are then in view of l|12|) and (|22[) 

$ (3) = ^— + me~ kz cos(kx) 

r\ 1 + 7] 

$W = 1 + ( v + e kz + v ~ e -kz\ cos ( fca; ) (24) 
7] 

— - + wie kz cos(fcx) . 

By using the linearized version of the boundary con- 
I 



ditions (|20|l and i|21|> we can express the amplitudes 
Mi, vf, , and Wx in terms of Ax and B\. This then 
allows us to expand the energy H19[) up to second order 
in Ai and B\ with the result 



£?( A a , B X ) (0, 0) + - 1 ^— - 2 H^k ' Q _ 2dk 1 + a fc^ 1 



P3 ?T e 

1/ ne~ 2dk — l \ e~ dk 

+ 2H! xt k ^ 2e _ 2dk _ 1 + k 2 j B\ + H^k (r, - 1) 7?2e _ 2dfe _ 1 ^iBi- (25) 

The energy has clearly a stationary point at Ax = Bi = 0. It is stable as long as the Hessian 

'I (p± - 2 Hi.t y;"!-, 1 , + „ t 2 ) (, - 1) 




is positive definite. An instability is signaled by a van- 
ishing determinant of "H. 

We note that, of course, exactly the same condition 
results from the usual procedure of linear stability anal- 
ysis. In this case one investigates the dispersion rela- 
tion oj{k) of interface deformations of the form = 
A\ exp(i(kx — cut)) and = d + Bx exp(i(kx — tut)) re- 
sulting from a linearization of the equations of motion. 
An instability occurs if u> acquires a positive imaginary 
part. The linearized equation of motion corresponds to 
the quadratic approximation of the energy. An advan- 
tage of the energetic approach is that it applies equally 
well to inviscid and viscous fluids. On the other hand 
it does not give information on the linear growth rate of 
the unstable perturbation. 

If the layer thickness d tends to infinity it can be in- 
ferred from 126|) that the off-diagonal elements of TC tend 
to zero whereas the diagonal elements reduce to the well- 
known form of a usual Rosensweig instability on a in- 
finitely deep layer of ferrofluid 0] . As expected we hence 
find in this limit two uncoupled interfaces showing inde- 
pendent Rosensweig instabilities of the usual kind. The 
situation is depicted in Fig. [21 where we have shown the 
determinant of the Hessian matrix TL as function of the 
wavenumber k for the case in which the critical field of 
the two instabilities coincides but the respective critical 
wavenumbers do not. 

The situation changes if the layer thickness is reduced 
as shown in Fig. [3J Due to the interaction between the 
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FIG. 2: Determinant of the Hessian for an infinitely thick 
ferrofluid layer as function of the dimensionless wave number 
k. Parameters pi = 2, ps = 0.5, a — 0.5, rj = 1.5, are chosen 
such that the two independent interfaces get unstable at the 
same value of the magnetic field, (by definition H c = 1), but 
at different wave numbers fc'"' = 1 for the upper interface 
and kc =2 for the lower one. 



surface deformations mediated by the magnetic field the 
degeneracy observed in the case d = oo is lifted and the 
lower layer "slaves" the upper one to its critical wavenum- 
ber. At the same time the critical wavenumber is shifted 
somewhat, k c ^ 1, from its "pure" value of the decou- 
pled case. The same holds true for the critical magnetic 
field strength. Moreover, the two interface deflections ac- 
commodate to each other in an anti-phase fashion. This 



5 



det (U) 





-0.02 




FIG. 3: Same as Fig. [5] for a layer thickness d — 2. The 
coupling between the two surfaces now lifts the degeneracy 
characteristic of Fig. [5] giving rise to new critical values for 
the wave vector, k c = 0.96, and the magnetic field, H c = 0.98. 
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FIG. 5: Dimensionless critical wavenumber of the linear in- 
stability as function of the ratio a of interface tensions for 
layer thickness d = oo (long dashed), d — 2 (dashed), d = 1 
(dotted), d = 0.5 (solid). Other parameters as in figB] 



manifests itself in different signs of A\ and B\ building 
the components of the eigenvector corresponding to the 
zero eigenvalue of TL. This anti-phase orientation was to 
be expected intuitively since it allows the largest gain in 
magnetic energy (cf. Fig[I}. 

Which interface dominates which depends on the pa- 
rameter values of the system and accordingly a crossover 
can be observed when some parameter is changed. In 
Figs. 14161 we give some examples of such crossover phe- 
nomena when the ratio a between the two interface ten- 
sions is changed. Fig. 0] displays the relative amplitude 
of the two surface deflections. The figure clearly indi- 
cates that for small values of a, i.e. when a^' <C <r^ d \ 
the lower instability dominates, A\ 3> B\, whereas with 
increasing a the amplitude of the lower interface deflec- 
tion decreases and the coupled unstable modes get more 
and more dominated from the upper interface. Similarly 
Figs. and show the crossover of the critical wavenum- 
ber k c and the critical field H c , respectively, when a is 
varied. In all cases the crossover gets sharper with in- 
creasing depth d of the fcrrofluid layer as expected. 




a 

FIG. 4: Relative amplitudes of the unstable modes as function 
of the ratio a of the interface tensions as defined in 11511 . The 
layer thickness is d=2 (dashed), d=l (dotted), and d=0.5 
(solid). The other parameter values are r\ = 0.66, pi = 1.2 
and p3 — 0.85. 

Although the linear analysis already reveals some as- 



1.05 



0.95 





upper interface 























0.4 0.5 0.6 0.7 0.8 0.9 1 

a 



FIG. 6: Dimensionless critical magnetic field H c of the linear 
instability as function of the ratio a between the interface ten- 
sions for layer thickness d = oo (long dashed), d = 2 (dashed), 
d = 1 (dotted), and d = 0.5 (solid). Other parameters as in 

figm 

pects of the interplay between the two Rosensweig in- 
stabilities it is not able to yield information about the 
static pattern of interface deflections that will eventually 
emerge. In order to address this problem we need to 
extend our analysis to include non-linear terms able to 
saturate the exponential growth predicted by the linear 
stability theory. This is the subject of the next section. 



IV. WEAKLY NONLINEAR ANALYSIS 

In our analysis of the energy of the surface deflections 
the instability of the flat surface corresponds to the mini- 
mum of the energy at A\ = B\ = turning into a saddle 
point. Within the quadratic approximation the energy 
hence decreases down to — oo with increasing amplitudes 
A\ and B\. In reality, however, already for moderate 
values of A\ and B\ higher order terms in the expan- 
sion of the energy have to be included which cure this 
divergence. As a result the energy again increases with 
increasing amplitudes A\ and B\ and correspondingly a 
new minimum forms describing the new stationary sur- 
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face profiles ( (a \x) and (( d \x). 

We assume that the susceptibility of the ferrofluid is 
sufficiently small such that an expansion of the energy 
(|19|l up to fourth order in the amplitudes of the surface 
deflection is sufficient to find the new stationary state. 
Such an expansion is equivalent to the derivation of a 
third order amplitude equation for the unstable mode 
|l5j . In order to obtain a consistent expansion the Fourier 
expansions ilL'Mli and (|24|1 have to be extended according 
to 



and 



$(3) 



C (d) (x) = d + J2 B n cos (nkx) 

n=l 

2 

f (°) ( X ) = A n cos (nkx) 



— h ^ u n c nkz cos(nkx) 

7] 1+77 



(27) 



TABLE I: Magnetic fluid parameters used in Figs. 17191 



Experimental parameters 


Dimensionless values 


= 1.69 g/cm 3 
p { > = 1.12 g/cm 
p (3) = 0.0013 g/cm 3 
cr (0) = 16.6 mN/m 
cr (d) = 25.9 mN/m 
d — 1.54 mm 
X = 0.8 


pi = l.5l 
pz = 0.001 

a = 0.64 

d = l.O 
r) = 0.29 



$ (2) = - + J2 ( v nZ nkz + v~c~ nkz ) cos(nkx) (28) 



n=l 
2 



= - + V w n e nfez cos(nfca;). 
»7 ±1 



FIG. 7: Stationary pattern of coupled surface deflections that 
evolve after the instability of the state with flat interfaces. 
The ferrofluid layer is shown in gray. Parameters are given 
in table Q the value of the external magnetic field is H ex t = 
l.OOOlf/c- The figure uses dimensionless units. 



To explicitly perform the minimization of the free en- 
ergya variant of the computer algebra code documented 
in is used. Fixing the desired order of the expansion 
(four in our case) this program selects in a first step those 
amplitude combinations which are compatible with the 
translational invariance of the problem in x-dircction. In 
our case only 17 of the originally 70 terms remain af- 
ter this procedure. In a second step the ansatzes (|27f) 
and (|28|l are used in the boundary conditions l|20() and 
(|21|l and the coefficients u n ,v£, v~,w n are determined as 
polynomials in the A n and B n . After this the energy 1)19(1 
can be expanded up to fourth order in A\ and B\ and 
up to second order in Ai and Bi- Several of the remain- 
ing terms disappear after the integration over x implicit 
in the horizontal average in (|19|) . Minimizing the result- 
ing expression in A2 and B2 we find that both are of 
order A\ , B\ which proofs the consistency of our expan- 
sion a-posteriori. Finally the free energy is minimized in 
the amplitudes A\ and B\ of the main modes. The final 
expressions arc explicit but too long to be displayed. 

For rf = oowe again reproduce the results obtained for 
the standard Rosensweig instability on a layer of infinite 
depth 0, ■ For d < 00 the two interfaces couple and 
the two surface deflections arrange in a stable, anti-phase 
pattern. 

To elucidate this final structure in detail we consider 
the case of experimentally realistic parameters collected 
in tabled From the linear stability analysis we find for 
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FIG. 8: Same as figE|for H cxt = 1.0Lff c . 
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FIG. 9: Same as %[7|for H cxt = 1.04iJ c . 
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the dimcnsionlcss wave number k c = 0.84 and for the cor- 
responding dimensionless field H c = 0.75. The station- 
ary interface profiles resulting from the weakly non-linear 
analysis are displayed in Figs. 17191 As can be seen the 
lower interface is the dominating one. For slightly over- 
critical magnetic field the lower interface already shows 
an array of developed Roscnsweig ridges whereas the up- 
per one is just gently curved by the inhomogencous mag- 
netic pressure resulting from the field modulation in- 
duced by the lower interface (Fig. 0. With increasing 
field both deformations grow (Fig. (SJ. A particular in- 
teresting case is shown in Fig. where the amplitudes of 
the surface deflections have increased to such an extent 
that the two interfaces touch each other. Correspond- 
ingly the ferrofluid layer stays no longer connected but 
disintegrates. In our two-dimensional (x, z) model this 
gives rise to the formation of parallel slices. In a more 
realistic three-dimensional setting, including surface vari- 
ations in y-direction as well, the layer would evolve into 
a regular array of disconnected islands. A similar phe- 
nomenon occurs in the usual Rosensweig instability on 
very shallow layers of ferrofluid Q. 

In Fig. 1101 we have shown the maximum and the min- 
imum layer thickness as function of the magnetic field. 
The formation of islands occurs when the lower branch 
intersects with the horizontal axis. Note that, at least 
for the parameters of table [IJ this happens already for a 
field exceeding the critical one by only 4%. 




0.735 0.74 0.745 0.75 0.755 0.76 0.765 0.77 
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FIG. 10: Maximum (full line) and minimum (dashed line) 
dimensionless thickness of the ferrofluid layer as function of 
the dimensionless external field strength for the parameters 
given in table Q The dots correspond to the situations dis- 
played in Figs. 17191 respectively. For the field at which the 
minimum distance between the interfaces shrinks to zero the 
layer disintegrates into an array of disconnected rolls. 

With the appearance of such an island structure our 
theoretical model breaks down. To study the future evo- 
lution of the structure when increasing the field still fur- 
ther it is more appropriate to start from a model of in- 
dependent ferrofluid drops [Tsj . 

In omitting higher orders of the expansion of the en- 
ergy in the amplitudes of surface deflection in our nonlin- 
ear analysis we have tacitly assumed that the fourth order 
terms are sufficient to saturate the linear instability, i.e. 
to make E[^°\ C (d) ] -> °° for ((C (0) ) 2 + (C (d) ) 2 ) oo. 



This, however, is correct only if the susceptibility \ of the 
ferrofluid is not too large and the thickness d of the layer 
is not too small. In Fig. ^] we have displayed the region 
in the d~x plane, in which our treatment is consistent. 



FIG. 11: Region of consistency of our nonlinear treatment of 
the pattern formation in the plane spannend by the dimen- 
sionless layer thickness d and the susceptibility x- Outside 
the shaded region the fourth order terms in the expansion of 
the energy are not sufficient to saturate the linear instability 
and higher order terms are needed to get finite results for the 
amplitudes of interface deflections when minimizing the en- 
ergy. The dashed line is the result of for a ferrofluid layer 
with rigid bottom. Parameter values are from table U] 



Finally Fig. 1 121 gives the phase diagram of the ferrofluid 
sandwich structure showing the transition lines from flat 
interfaces to anti-phase interface modulations and fur- 
ther to disconnected regions. It would be interesting to 
compare the location of these theoretical lines with ex- 
perimental results. 



I / 


II 


= ^^^^^ 1±± = 



0.68 0.7 0.72 0.74 0.76 0.78 0.8 

Hat 

FIG. 12: Phase diagram in the plane spanned by the dimen- 
sionless external field 7i7 cxt and layer thickness d for a fer- 
rofluid sandwich structure with the parameters given in table 
H] In region I both interfaces are flat, in region II, ridges oc- 
cur and in region III, the layer disintegrates into an array of 
disconnected rolls. The dashed line indicates that for larger 
values of the magnetic field higher order terms are necessary 
to accurately determine the location of the transition line. In 
the shaded region the fourth order terms in the energy are 
not sufficient to saturate the linear instability and higher or- 
der terms are mandatory, cf. Fig. 1111 
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V. CONCLUSION 

In the present paper we have investigated the linear 
and weakly nonlinear theory of two coupled Roscnsweig 
instabilities in a fcrrofluid sandwich structure. To this 
end an approximate expression for the energy of the sys- 
tem was minimized in the deflection amplitudes of the 
two interfaces between magnetic and non-magnetic liq- 
uids. The approximate expression for the free energy 
was obtained from a fourth order perturbative expansion 
in these interface deflections. At the onset of instability 
the two individual Rosensweig instabilities compete and 
depending on the concrete values of the parameters one is 
able to "slave" the other one to its unstable wavenumbcr. 
As a result, a stable antiphase pattern of two interacting 
modulated interfaces arises. For sufficiently thin layers 
and sufficiently large magnetic fields the two curved in- 
terfaces may touch each other which brings about the 
disintegration of the layer and gives rise to disconnected 



rolls or islands. Using realistic parameter values we gave 
estimates of the required layer thicknesses and magnetic 
fields necessary to observe this phenomenon in an exper- 
iment. Being perturbative in nature our theoretical anal- 
ysis has a limited range of validity which we quantified by 
estimating the contributions of higher order terms. It is 
possible though tedious to push the expansion to higher 
orders in a systematic way. 
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